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Abstract 

Geometrical random multiplicative cascade processes are often used to model positive- 
valued multifractal fields such as for example the energy dissipation field of fully 
developed turbulence. A dynamical generalisation of these models is proposed, which 
describes the continuous and homogeneous stochastic evolution of the field in one 
space and one time dimension. Two-point correlation functions are calculated. 
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Whenever strong anomalous, intermittent fluctuations, long-range correla- 
tions, multi-scale structuring and selfsimilarity go hand in hand, the label 
'multifractality' is attached to the underlying process. Although we have fully 
developed turbulence of fluid mechanics [1] in mind, such processes occur in 
various diverse flelds such as formation of cloud and rain flelds of geophysics 
[2], internet traffic of signal and communication engineering [3], and stock re- 
turns of econometrics [4], just to name a few. The maybe most widely used 
textbook [5] visualisation for multifractality is represented by a multiplicative 
cascade process. It reproduces the above mentioned properties by introducing 
a scale-independent cascade generator, which produces a nested hierarchy of 
scales and redistributes the local mass multiplicatively. 

In fully developed turbulence, random multiplicative cascade models (RMCM) 
are often employed to model the energy cascade, describing the energy flux 
through inertial range scales. Due to their multiplicative nature, it is straight- 
forward for them to reproduce multifractal scaling exponents associated with 
the energy dissipation field [6], the latter representing the intermittency cor- 
rections to the K41-thcory [1]. Although the justification of RMCMs in terms 
of the Navier-Stokes equation is far from being clear, these phenomenological 
models appear to contain more truth than originally anticipated, as recent 
investigations on multiplier distributions [7] and scale correlations [8] have 
revealed. - On the other side RMCMs are purely geometrical in nature and 
are not capable of describing causal dynamical effects of the turbulent en- 
ergy cascade like backscattering of energy fiux from small to larger scales. A 
generalisation in this direction is clearly called for. 

In this Letter we present such a generalisation and construct a dynamical 
RMCM in 1 + 1 space-time dimensions, respecting causality. Multifractal scal- 
ing is recovered for n-point correlation functions. 

In a discrete RMCM in 1-1-0 dimensions, the amphtude 



of the positive- valued energy-dissipation field, resolved at the dissipation scale 
Tj, is equal to the product of independently and identically distributed ran- 
dom weights q{lj), which come with (q) — 1 and a nested hierarchy of scales 
V = l.j < Ij = L/\^ < Iq — L with Q<j<J. The integral length L and the 
dissipation length t] represent the largest and smallest length scale involved, 
whereas A>1 is the discrete scale step. 

It is the second step of (1), which we now exploit to propose a generalisation 
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of the cascade field, now defined in 1+1 space-time dimensions: 

(oo oo "\ 

J dt' J dx'f{x -x',t- t')-f{x', t') \ . (2) 
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Due to the point-resolution the function 7(x, t) ~ Sa{{dxdt)^ ^^a, —\, ji) is 
taken as a stable white-noise field with index 0<«<2 [9]; for the special case of 
a=2 it corresponds to a non-centred Gaussian white-noise field. Its character- 
istic function is (exp{n7}) = exp{— ((T"n")/(cos(7rQ;/2))-|-//n}, which fixes the 
parameter = cr"/ cos(7rQ;/2) in order to fulfill the expectation (exp{7}) = 1. 
Imposing causality, i.e. the field amplitude t) is only influenced by the 
past and not by the future, requires /(x, t) = for i < 0. This is satisfied with 
the following symmetric index function 



1 {^<t<T, -g{t)<x<g{t)) 
(otherwise) , 



where the window function g{t) introduces a correlation time T and a corre- 
lation length L with g{Q) a; and g{T) = L/2] see Fig. 1 for an illustration. 
The exponent of the construction (2) can be thought of as a moving average 
over the stable white-noise field. 

According to (3), the time integration in (2) goes from —T<t'<S), where x — 
t — Q has been set for simplicity. Upon assigning 2g{—tj)—lj the hierarchy of 
length scales Ij is translated into a hierarchy of time scales, so that (2) then 
factorises into contributions 



q{lj) — exp < 



' *i 9i-t') 

J dt' J dx'j{x',t') 

tj-1 -g{-t') 



(4) 



see again Fig. 1. In order to interpret this as a random multiplicative weight, 
the probability distribution density of q{lj) needs to be independent of scale. 
Since the j{x', t') are i.i.d., the integration domain of (4) then has to be inde- 
pendent of the scale index j. This fixes the window function within —T<t'< — 

a(-t') - ^^'^^ (5) 

+ „ {T-tr,) 



which also satisfies the boundary conditions g{T)=L/2 and g{tn)=r]/2. So far, 
no specification of g{—t') for — < i' < has been given; since tj^<^T should 
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hold on physical grounds, the simplest choice would be tj, — 0. 

The construction proposed in Eqs. (2)-(5) guarantees that the one-point statis- 
tics of the dynamical RMCM is identical to its geometrical counterpart. In 
order to qualify for a complete dynamical generalisation, it is not only the 
one-point statistics, but the n-point statistics in general, which should match. 
Hence, we now consider the equal-time two-point correlator 



(6) 



The correlation between the two points with distance 1<L stems from the 
overlap region of the two index functions f{xi, t) and f{x2=Xi + l,t)] see Fig. 
2a. This explains the second step of (6), where 



, -9-l(Z/2) 9{-t') , 

D{l>ri)=expl J dt' J dx'-f{x',t')\ (7) 

i-9{-t') 



represents the contribution from the overlap region. The contributions from 
the non-overlapping regions are statistically independent, hence factorise and 
cancel. Introducing the spatio-temporal overlap volume 



-g-l(//2) g(-t') 

V{l>T])= J dt' J dx' 



(L-n) \U (L-n) 

and employing basic properties of stable distributions [9], the expectation of 
the n-th power of expression (7) can be transformed into 

{D{ir)^exp(^V{l){n-n'^)] . (9) 

\COS 2 / 



Defining the multifractal scaling exponents r(n) = r(2)(n— n")/(2— 2") with 
r(2) = (a"/cos^)(2-2")r/L(T-t^)/(L-?7) as well asr[ni,n2] =r(ni+n2)- 
Ti'iT'i) —'^{'^2), insertion of (9) into (6) leads to the final expression for the two- 
point correlator: 
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It reveals niTiltiscaling behaviour for ri<l<^L, which shows that the equal-time 
two-point statistics of the dynamical RMCM in 1-|-1 dimensions is completely 
analogous to the findings of the geometrical RMCM in 1-1-0 dimensions. 

In view of this analogy, the large-scale deviation appearing in the last expres- 
sion is reminiscent of similar recent findings on scaling functions within the 
geometrical RMCM [10]. With the introduction of a suitably tuned large-scale 
fiuctuation this large-scale deviation can be cancelled. In our present context 
this implies a small extension of the window function (5) beyond the time 
interval -T<t'<0: g{T<-t'<T+AT) = L/2; consult again Fig. 1. It is then 
straightforward to show that for AT = (T — t^)ri/{L — 1]), which is much 
smaller than the correlation time T, the two-point correlator exactly becomes 
Cn„nAl) = {L/iy^^^^^^^ for all 7]<1<L. 

Besides equal-time two-point correlations temporal two-point correlations are 
also of interest. These correlations arise due to a modified spatio-temporal 
overlap volume as illustrated in Fig. 2b. A straightforward calculation, anal- 
ogous to (7)-(9), yields 



t n ( i\\-r\nxM 
T-ix T[ni,n2] 

l) ■• 

for simplicity, the parameter of the window function (5) has been set tr, — 0. 
The last step, only holding for 0<^t<T and r)<^L, reveals that the temporal 
two-point correlator comes with scaling exponents identical to those of the 
equal-time two-point correlators. 

Within this dynamical RMCM investigations on several other observables are 
nearly as straightforward [11]. Here we mention only spatio-temporal (n>2)- 
point correlations in connection with fusion rules, and moments and multiplier 
distributions based on coarse-grained field amphtudes. All of these observables 
are found to be in good agreement with experimental results derived from the 
surrogate energy dissipation field of fully developed turbulence. 

The dynamical RMCM, presented here, is a generalisation of the geometrical 
RMCM. By construction it is continuous and homogeneous, does not make 
use of a discrete hierarchy of scales and stochastically evolves a positive- valued 
field in one space and one time dimension. Several generalisations of this new 
model immediately come to mind: stochastic evolution in n+1 dimensions 
with the optional inclusion of spatial anisotropy, discretisation of space-time 
into smallest cells to model dissipation and deviation from log-stability, and 
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dynamical RMCM for vector fields to model the turbulent velocity field. Work 
in these directions is in progress. 
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Figure 1: Illustration of the causal index function (3) used to construct the 
positive- valued multifractal field e{x,t). 
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